Distorted wave impulse approximation analysis for spin observables in nucleon 
quasi-elastic scattering and enhancement of the spin-longitudinal response 
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o : 

' We present a formalism of distorted wave impulse approximation (DWIA) for analyzing spin observables in 

^ nucleon inelastic and charge exchange reactions leading to the continuum. It utilizes response functions calculated 

by the continuum random phase approximation (RPA), which include the effective mass, the spreading widths and 
^ ' the A degrees of freedom. The Fermi motion is treated by the optimal factorization, and the non-locality of the 

O . nucleon- nucleon t-matrix by an averaged reaction plane approximation. By using the formalism we calculated the 

spin-longitudinal and the spin-transverse cross sections, IDg and IDp, of ^^C, ^°Ca (p, n) at 494 and 346 MeV. 
The calculation reasonably reproduced the observed IDq, which is consistent with the predicted enhancement of 
the spin-longitudinal response function i?L. However, the observed IDp is much larger than the calculated one, 
which was consistent with neither the predicted quenching nor the spin-transverse response function Rt obtained 
by the (e, e') scattering. The Landau- Migdal parameter gjvA for tbe A^'A transition interaction and the effective 
^ ■ mass at the nuclear center m* (r = 0) are treated as adjustable parameters. The present analysis indicates that the 

' smaller g^^i^ 0.3) and m*(0) « 0.7m are preferable. We also investigate the validity of the plane wave impulse 

y—^ , approximation (PWIA) with the effective nucleon number approximation for the absorption, by means of which 

' i?L and Rt have conventionally been extracted. 

00 ; 

O . PACS numbers; 24.70. -f-s, 24.10.-i, 25.40. Kv, 21.60.Jz 

o ■ 

■ 

^ ; I. INTRODUCTION 
^ '. 

^ , The study of nuclei by means of intermediate energy (~100 MeV- ~1 GeV) proton beams has been very active 

■ since the end of the 1970's, and has led to various advances such as the discovery of the Gamov- Teller (GT) giant 

] resonance and finding the quenching of their strength. In the 1980's great progress was made in experimental facilities, 

• ^ which can afford to accelerate polarized proton beams with intermediate energies and to measure the polarization of 

K% scattered protons and neutrons and so on. This made it possible to carry out complete measurements of {p,p') and 

^ . (p, n) scattering, namely, measurement of the polarization P, the analyzing power Ay and the polarization transfer 

- - ' coefficients Dij . Theoretical investigations related to these experiments, such as search for the origin of the quenching 
of the GT strength and studies of the precursor phenomena of the pion condensation, have also been pursued vigorously 

m- . o 

In the course of these activities, a very interesting prediction was presented by Alberico et al. H at the beginning 
of the 1980's. They claimed that in the quasi-elastic region with fairly large momentum transfer 1-3 fm^^) the 
isovector spin-longitudinal response function Ri^{q,u!) is enhanced and softened while the isovector spin-transverse 
response function i?T('?,w) is quenched and hardened, where w is the transferred energy. The enhancement of i?L 
is attributed to the collectivity induced by the one-pion exchange interaction and thereby understood as one of the 
precursor phenomena of pion condensation. The quenching of Rt, on the other hand, is due to the combined effect 
of the repulsive short-range correlation and the one rho-meson exchange interaction. 

A great deal of experimental work has been done in order to explore this prediction. Measurements of the po- 
larization transfer coefficients Dij were carried out at LAMPF for {%>,%>') at an incident energy of 494 MeV, from 
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which the ratios i?L I Rt were extracted ■ Surprisingly, the ratios were less than or equal to unity, which seriously 
contradicted the prediction i?L 3> i?T- However, the scattering {p,p') mixes the isoscalar and isovector contributions, 
and thus the estimation of the ratios was not conclusive. 

Later complete measurements of the (p, n) reaction, which focused exclusively on the isovector contribution, were 
carried out at LAMPF and at RCNP ||]. The ratios obtained were still less than unity. From these results, it 
was concluded |lOyil| ] that there is no enhancement of i?L, namely no collective enhancement of the pionic modes, 
which was interpreted as evidence against the pion excess in the nucleus. With help of sum rules, Koltun Jl^ analyzed 
the data by means of a correlated nuclear theory, in which the two-nucleon correlation is dominating, and concluded 
that there is no disagreement between the data and the model for i?L, but Rt is not explained. He claimed that there 
is no collective enhancement of the pionic modes. However, we believe that there are many questions to be solved 
before such conclusions are reached. 

The experimental evaluation of the ratio was based on the plane wave impulse approximation (PWIA) with effective 
nucleon number approximation for absorption where the number of participant nuclcons is estimated by a simple 
eikonal approximation. Evidently we should take into account the effects of distortion more accurately. Since this 
causes the nuclear responses to be of surface nature, the semi-infinite slab model |ll3| , |l4| , the surface response model 
[ p5| , p^ and others, were introduced. The distorted wave Born approximation (DWBA) analysis of the continuum 
spectra of the intermediate energy (p, n) reaction was carried out by Izumoto et al. [|l7| where the response functions 
were calculated by the local density approximation. The calculation in the distorted wave impulse approximation 
(DWIA) for the spin- longitudinal and -transverse cross sections, IDq and /-Dp, in the continuum was first performed 
by Ichimura et al. p8[ | for (p^jf). We utilized the response functions calculated by the continuum random phase 
approximation (continuum RPA) with orthogonality condition including the A degrees of freedom. We there 
employed the optimal factorization prescription pO|-p^ to deal with the Fermi motion of the struck nucleon, restricted 
the spin dependence of the nucleon- nucleon (iViV) t-matrix to the terms with ((Tq • q)(cri • q), and (ctq x q) • (<ti x q) 
and neglected the interference between the different spin-dependent terms of the NN t-matrix. De Pace |Q treated 
the absorption by the Glauber approximation and calculated the spin-longitudinal and -transverse cross section up 
to the two-step processes. He also used the same assumption for the NN t-matrix. Recently, Kim et al. also 
developed a DWIA formalism for calculating the spin observables in a form of inhomogeneous coupled-channel integral 
equations in the Tamm-Dancoff approximation. They also used the same assumption for the NN interaction and 
further neglected the spin orbit force in the optical potential. Noting this situation, we definitely need a more reliable 
method for analyzing the spin observables as well as the inclusive cross sections in the continuum. 

In this paper, we develop a DWIA formalism for nucleon-nucleus {NA) scattering at the intermediate energy leading 
to the continuum, with the response functions non-diagonal with respect to the momentum transfer as well as the spin 
directions. This method does not require the above restriction on the NN i-matrix and can handle the interference 
between the different spin-dependent terms. Thus it is much more reliable for calculating the spin observables such 
as P, Ay and Dij . 

The original prediction of Alberico et al. ^ was based on the theoretical framework of (1) the Fermi gas model, 
(2) RPA, or more precisely the ring approximation, with the one-pion -I- one-rho-meson exchange interactions -I- the 
contact interaction specified by the Landau-Migdal parameters, (/"s {ir + p + g' model), and (3) the universality ansatz, 
9nn ~ 9na ~ .9aa- ^ number of improvements have been made to various aspects of the method, such as use of the 
continuum RPA to treat the nuclear size effects [^^, and removal of the universality ansatz 27|. In the previous 
continuum RPA calculations, the mean field was assumed to be local and the spreading width of the particle was 
treated by the complex potential, but that of the hole was neglected. In this paper we take account of the non-locality 
of the mean field by the radial-dependent effective mass, m*(r), and the spreading width of the hole by the complex 
binding energy. 

Special importance has been put on the choice of the Landau-Migdal parameters. The position of the GT resonance 
gives severe restriction on the value of g^jv (~ 0-^ ~ 0.7). Effects of A-isobar are very sensitive to the value of (7^^ 
|g6|,^, which is also crucial for the cause of the quenching of the GT strength |^^. In the present analysis of the 
quasi-elastic (p, n) reaction, we treat g'pf^ and the effective mass at the center of the nucleus m* (r = 0) as the 
adjustable parameters, and thus try to evaluate the value of g'j^^. This is one of the main aims of the present paper. 

This paper is organized as follows. In Sec. |l| we summarize the general formulas for the cross sections and the spin 



observables of the {N, N') reactions, with special attention to the polarized cross sections IDi. In Sec. IH, the PWIA 
formalism with the optimal factorization and the effective nucleon number approximation is briefly reviewed because 
it has been used to extract i?L and Rt up to now. 

In Sec. 1^ we present the DWIA formalism for the intermediate energy (N, N') reactions to the continuum in 
cooperation with the response functions, again utilizing the optimal factorization. Since we treat the finite nucleus, 
calculations are carried out in the angular momentum representation, details of which are presented in Sec. 0. We 
describe a way of calculating the response functions non-diagonal in the coordinate space as well as the spin space. 
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which involves the radial-dependent effective mass and the spreading widths of both the particle and the hole in 



Sec. VI 



In Sec. VII, we perform numerical analysis for the reactions ^^C, ^'^Ca{p,n) at 494 MeV and at 346 MeV. We 
see that the spin-longitudinal cross sections are reasonably well reproduced. This is consistent with the predicted 
enhancement of i?L though the ratio Ri^/Rt is less than unity. However, the spin-transverse cross sections are very 
much underestimated. The contradiction regarding the ratio seems to come from the large difference between the 
experime nts an d the theories in I Dp. This confirms the reported conclusion in the experimental papers 



In Sec. VIII we test the reliability of the conventional method for extracting the response functions by comparing 



the results of DWIA and PWIA with the effective nucleon number. Sec. IX is devoted to some other questions such 
as the effects of the spin-orbit force and the ambiguity of the optical potential. Sec. |^ consists of a summary and 
conclusion. 

II. GENERAL FORMULAS FOR CROSS SECTION AND SPIN OBSERVABLES 

We consider the nucleon-nucleus (NA) inelastic and charge- exchange reactions to the continuum. First we sum- 
marize the general formulas for the double differential cross section and the spin observables of a polarized nucleon 
scattering off a nucleus, 

N + A-^ N' + B. (2.1) 

In the NA center of mass (cm.) system, we denote the scattering angle by Ocm, and the momenta of the incident 
and outgoing nucleons by k; and k/, respectively. The momentum transfer to the scattering nucleon is given by 

q = k/-k„ (2.2) 

and the energy transfer to the target is written as 

iJcm ^EN{k^)^EN'{kf), (2.3) 



where Ea{k) = + rrv^ is the energy of a particle a with mass TOq,. We use the unit system fi ^ \ and c = 1 
throughout this paper. 



As the coordinate system, we use either 

k 

M " |k. xk^l 



^ = ^^> y = ^^^. x^yxz (2.4) 



q » X kf , 

These are called the x, y, z and g, n,p directions, respectively. 

The unpolarized double differential cross section is expressed by the N A t-matrix T as 

I {9 cm, l^cm) 



2(2Ja + 1) (27r) 

x5iLJcn.-{EB{kf)~EAih))), (2.6) 

where rrig.ijnsj) is the spin projection of the incident (outgoing) nucleon and J a is the target spin, and /ij and /xy are 
the relativistic reduced energies 

EN{k{)EA{h) EN'{kf)EBikf) 
Mi = 7= , /^/ = -7= — — (2-7) 

^Jsna ^Jsna 

with SNA = {EN{ki) + EA{ki))^ . The wave functions 5'„o and are intrinsic states of the target A and of the 
residual nucleus B, respectively. They are governed by the intrinsic Hamiltonian Ha of the A-body system as 
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Ha^u = £n^n (2.8) 

where Ha includes the mass terms and thus means the invariant mass of the ^-body systems. Note that the final 
states are mostly unbound. In the summation S^^, uq runs only over the degenerate ground states of the target 
A. 

To extract information of intrinsic states separately, we introduce the intrinsic energy transfer lu and rewrite the 



cross section of Eq. (2.6) as 

/(0c„.,^cm) = 2i2jl+l) ^' ' '^^■^^ 

with 

Tr' [TTt] = ^^'(vI/„|T|*„„)(vl/„jrt|vI/„)5 (^ - (f„ - f„J) , (2.10) 

n no 

where Tr means the trace of the spin states of the incident and exit nucleons. The kinematical factor K is given by 

_ MiM/ kf duj ^ n.fif kf ^/JWa , . 

(27r)2 fc, (27r)2 k, ^ ' ' 



using the relation = ||3l|l. We note that the present K is equal to 2(2Ja + 1)CX of Ref. g , 2(2 Ja+I)^^^^ 

of Ref. |§ and of Ref. 



We represent the Pauli spin operators in the i-direction of the fc-th nucleon by aki- To unify the notation we 
also introduce ako{= 1/c) for the unit spin matrix of the fc-th nucleon. The nucleon number fc = denotes the 
incident or exit nucleon. Then the polarization, the analyzing power and the polarization transfer coefhcients are 
given respectively by 

p ^ TTTT'[TT^aoy] ^ ^ TrTr^jTap^Tt] 
^ TrTr'[Trt] ' ^ TrTr'[rrt] ' 



TrTyjTVoiTVoi] 
TrTr'[rrt] 



A. = (2.12) 



where «, j — x,y,z or q, 

The DWIA calculation is usually carried out in the [x, y, z] frame, while the [q, n, p] frame is sometimes more 
convenient for the theoretical analysis. The relation between Dij — q, n,p) and Dij = x, y, z) is given by 



Dyyi 



( 


Dp,] 


-i 




D,J 





^ zz ^ zx 

Smt/n COS t/n I \ O ry D^r J \ Sill t^n COS ( 



(2.13) 



where dp is the angle between p and z. 

To extract nuclear responses, Bleszynski et al. ||2^ decomposed the NA i-matrix in the [q, n, p] frame as 

T = Tocrao +Tn(TQn +Tq(7og +Tpaop (2-14) 

and introduced the polarized cross sections IDi, which extract Ti exclusively as 

I Do = ^[1 + + D,, + Dpp] = ^j^Tr' [ToT^], 

ID,, = i[l + Dr,„- - Dpp] = _i^Tr'[r„T^], 

ID, = ^[1 - D„„ + D,, - Dpp] = -J^Tr'[r,Tt], 

4'' "^'^ ' 2Ja 



I Dp = - [1 - Am - A? + Dpp] = , , Tr' [TpTj] . (2.15) 
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The unpolarized cross section / is expressed as 

I = IDa + IDn + IDq + I Dp. 



(2.16) 



In the N A laboratory system we denote the angle, the momenta and the energy transfer corresponding to 0cm, k^, 
and Wcm by 6'iab, Ki,K/ and wiab, respectively. The unit vectors 



N = N' = 



K, X K 



/ 



|K, xK^r 



L = K„ S = NxL, 



L' = K/, S' = N'xL' 



are usually used to specify the directions, and are denoted by iV, L, 5*, N' , L' and S", respectively. 
In this system, the unpolarized cross section is given by 



'lab^C'lab 



^^lab) — 



dfliahdcjiah kf 



due to the relation l32| 



sm( 



dillabC?l^lab sin (?iab 



and the observed polarization transfer coefficients Dij{i — S,N,L;j = S' , N' , L') are obtained from the 
Dij{i,j = x,y,z) as 



Dnn' = D 

f Dll' Dls 
\DsL' Dss 



yy 

^ zz ^ zx 



coscki — smai 
sinai cosai 



(2.17) 

(2.18) 

(2.19) 
calculated 

(2.20) 



with ai — 0iab + r ^, wh er e fl i s the r elatiy istic spin rotation angle [ pO| . 

By use of Eqs. ( 2.13 ), ( ^.15| ) and ( 2.20 ), Di^s are obtained from -Dy 's in the NA laboratory frame as 



- + Dnn + {Dss' + Dll') cos ai 

+ (DsL' - Dls') sinai], 
Dn = T fl + Dnn - {Dss' + Dll') cos ai 

- {Dsu - Dls') sinai], 
NN + {Dss' - Dll') cos a2 

- {Dls' + DsL') sin 02], 

NN - [DsS' - DlL') C0SQ;2 

+ {Dls' + DsL') sin 0^2] 



4 
1 

D, = \[l-D 
D„ = -\l-D 



(2.21) 



with a2 = 29p — ai. 



III. PWIA FORMALISM 

Before presenting the DWIA formalism, we briefly review the PWIA formalism ||3^] , which has conventionally been 
used in the analysis of quasi-elastic scattering to extract the spin response functions pHsl. To avoid confusion, we 



suppress the isospin and A degrees of freedom until Sec. HID, where we discuss these degrees of freedom in detail 
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A. PWIA i-matrix 



The PWIA i-matrix in the NA cm. system is written as 

A 



Tmkf,k,) = (M/„|5^tfc(k/,k,)|*„o), (3.1) 

k=l 

where ifc(ky,ki) is the NN scattering t-matrix between the incident nucleon and the fc-th nucleon in the nuclei. 

To avoid the difficuhy of Fermi momentum integration, the optimal factorization approximation |^-^,^ is often 
used, where t/c(kj,ki) is replaced by the NN t-matrix in the optimal frame, t5?(kj,k,;), which is written as 

*fc(k/,k,) = (k/,p'|t|k,,p) ="^tlf^Ckf,k,)aoaCrkb 

ab 
n 

+ D'la^pCfkq + D^cjQqakp + E^aoqUkq + F^aopfTkp- (3.2) 
The optimal momenta p and p' of the struck nucleon in the nucleus are given by 

p-(^-^/)q-^^^, p' = p-q, (3.3) 

respectively, and the parameter rj is determined by the on-shcU condition £'Ar(k,;) + En{p) = _Ejv(k/) + i?jv(p')- 
This NN t-niatrix in the optimal frame is obtained from the observed NN t-matrix in the cm. frame, 



t™(K',K) = {K\~K'\t\K,-K) = ^C,'J,(/«',/«)aoa'(Tlb' 



a',b' 



= A'ctooCTio + B'cFQn^ain^ + C"(cro„„Crio + CTooCTlnJ + E'cFQq^aiq^ + F'aQp^aip^, (3.4) 

where k {k') is the initial (final) relative momentum in the NN cm. scattering, which is determined from k^ (kj) 
and p (p'). Here the coordinate system is determined by 

n, Pc = Qc X ric (3.5) 



Iqd \ KX K' 

with qc = k' k. The relation between i''(k/, k;) and t™{n' , k) is given by 

i''(k/,k,) = J,(k^,k,)i?fpi„r°(/^',/^)i?^pi„. (3.6) 

The MoUer factor is given by 

J.(k/, k,) = (3.7) 

where we neglect the mass difference between the proton and the neutron. The relativistic spin rotations R^p^^ and 
i?gpijj are given in Eq. (3.29) of Ref. and the relation between and t™, is explicitly given in Eq. (3.34) of 
Ref. |0). 

The approximation q w qc holds [pO| for the reactions we are interested in. Noting that fie = n, we can identify the 



coordinate system [q, fi, p] with [qc,nc,Pc]. We also showed [|0| that the approximation RJ^^^^ ~ 1 holds. Therefore 



4, 11, pj wiuii i^nc , lie , fcj • vvc aiDu oiiuwcu | pu| uiidu uiic ajjjjiuT^iiiiauiuii 

we can safely write 

t:,(kf,k.o = ukf,k,)t:'^iK',K), (3.8) 

namely, 

Dl^Dll^ 0, E^-i = Jr,E\ F" = Jr,F'. (3.9) 



From Eqs. (BjJ) and (3.2), the PWIA t-matrix with the optimal factorization is given by 
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TPW(k/,k,J ^ ^^i^^(k/,k,)Fr°(q)^Oa (3.10) 



a b 

with the transition form factor 



Fr°(q) = (*„|^c7fe,e-'''--'=J rfe,] |vI/„„). (3.11) 



B. Response Functions 



From Eq. (2.15), the polarized cross sections /lab^'i in PWIA are expressed as 

iriDf^ - ifiabEC(k/,kOi?.a(g,c.)tI(k/,kO, (3.12) 



where Kui, — {Kf/kf)K (see Eq. ( ^.19 )) and Ri,a{q,uj) are the response functions 



Rba{q,u;) = ^^^i-^^^Vr"(q)Fr''*(q)5(t^- (f„ -f„o)), (3.13) 

n no 

which are determined solely by the nuclear intrinsic states. Note that Ri,a depends on only the magnitude q = |q|, 
because of the unpolarized target, i.e. the presence of the summation 

For comparison with the results of different nuclei, the normalized response functions 

Rba{q,i^) = ^Rba{q,uj), (3.14) 

whose spin diagonal parts satisfy the relation 

Raa{q,(^)duj ^ 1 for g ^ oo, (3.15) 



are more conve nien t. 

Noting Eq. (3^), I\'J^D^^ and -^lab^-D^^ are simply written, using the normalized spin- longitudinal and spin- 



transverse response functions i?L and Rt as 



ir^^D-^ = K,,^A\E^kf,k,rR^iq, 



IilD^p"" = Ki,^A\F^ikf,K)\'RTiq,uj) (3.16) 
where 

Rh{q,i^) = Rqq{q,Uj), 

R-riq.uj) = - {kpp{q,uj) + Rnn[q,i^)^ ^ Rpp{q,uj). (3.17) 

C. Effective nucleon number approximation for absorption 

In practical analysis we must include the effects of the absorption, which are usually treated by the effective nucleon 
number Agff estimated by the simple eikonal approximation as 

/•OO 

Acs = / 27r6n(6) cxp [— n(6)t7Ar7v] d6 (3.18) 
Jo 



with n{b) = dzpA{\^ + &^) where h is the impact parameter, ctattv is the total NN cross section in the nuclear 
medium, and pA is the nuclear density. 
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Then Jiab^^i's in PWIA with absorption, hahDi, are given with A^s by 

ah 

and especially 

/lab^g = iflabAeff|i;''(k/,k,)|'i?L(g,C^), 
/lab^p = iflabAeff|F''(k/,k,)pi?T(<Z,^). 

This prescription is called the effective nucleon number approximation. 



(3.19) 



(3.20) 



D. Isospin and A degrees of freedom 



Now let us take account of the isospin and the A isobar degrees of freedom. Following Ref. |27| with some 
modifications, we introduce the unified notation of spin and isospin operators, crj^^j and r^^ , of the fc-th particle [N or 
A) as 



and 



CTfea for a ~ N 
^ka = \ Sl^ for a = A 
Ska for a = A* 



for a = N 
(tI), fora = A 
Tk^ for a = A* 



(3.21) 



with K = s, 0, +, — where 



Tfcs — Ifc, Tk± 



V2 



, TkO — Tkz 



(3.22) 



(3.23) 



and Sfc and Tfc are the standard spin a nd i sospin transition operators from N to A. 

We extend the NN t-matrix of Eq. (3.2) to that which includes the isospin and the A^A channel in a restricted 
form as 



,r,_,>)(0) ,77(1)^ ,V(A) rp , ,»7(A*) rpf 

tk-t^ +tf, TQ-Tk+tf^ To-Lk+tk To- If, 

a=Af,A,A* K=s,0,± ab 



(3.24) 



Note that = (T^)]: = r^f'^'^ and so on from the definition ( |3.22| ). The PWIA t-matrix of Eq. (|3.10| ) now becomes 



TP„^(k,,k,) = EEECab(k/,k.)F:--(q)«, 

a K ab 

where the transition form factor of Eq. ( |3.11 ) is generalized as 



(3.25) 



fc=l 



Then the response function of Eq. ( 3.13 ) becomes 

n no 

where only i?os and i?sO remain finite due to the charge conservation. 
Now the polarized cross sections Jiab-Dj are expressed in PWIA as 



(3.26) 



(3.27) 
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rPW 
Mab 



><EECa(k/,kOi?L%a(9,'^)C(k/.kO- (3.28) 

a/3 afo 

The normalized response functions are usually introduced only for the diagonal parts in the isospin space as ^ 

1 



(3.29) 



with 



A/L = 2iV, A/'o = A4 = ^, A/V = 2Z, (3.30) 
where N and Z are the numbers of neutrons and protons, respectively. The normalization relation then becomes 

R^.Z{q,L^)du;^l for g ^ oo, (3.31) 



if the A component does not exist in the ground state. 

To use Eq. (3.28), we need information about t^^ as well as i**^, but to our knowledge it is not available for such 
a very off-shell region. Suggested by the n + p exchange model, we assume that 



,?)A _ /a ,rjN 
JN 



(3.32) 



for K = 0, +,— and i = q,n,p where /jv and /a are the ttNN and the TriVA coupling constants, respectively. Then 
we define the normalized isovector spin-longitudinal and -transverse response functions with the A contribution, i?K,L 
and Rk.,t, as 12^ 



-[l,K,L — J^K,,qq ' ^ I /jv ' ''''' 
jNN , j Ia \ f,AN 



JN ' 



N 



/a \ AAA 



(3.33) 



in place of Eq. ( |3.17| ) . 

For the {p,n) reaction, Eq. ( 3.20| ) are generalized as 



/lab^, = 4XiabA^cff|£^''(k/,k,)p7?„i(g,L^), 



(3.34) 



by using Eqs. (|3!2^)- (^!30|), (^3^ ) and {n\TQ^\p) {p\T^_\n) = 2. The effective neutron number Ncs is defined in a 
similar way as that of Eq. (3.18). For the {n,p) reaction, N^g should be replaced by the effective proton number Z^s 
and it!_ L(x) by i?+ L(T)- For the isovector part of the {p^p') or {n,n') scattering, 4iVoff should be replaced by A^s 
and -R-,L(T) by -Ro,l(t)- 



IV. DWIA FORMALISM 



We now present ou r DW IA formalism. To avoid unnecessary complexity, we suppress the isospin and the A degrees 
of freedom until Sec. VI D 



A. DWIA t-matrix in momentum representation 

The matrix elements of the DWIA i-matrix in the NA cm. system are written as 
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[?^°Y(k/,k.)]„^ - (k')*„ix:ifc(k',k)i*„„xi:ijk)) (4.1) 

^ ' fc=i 

in the momentum representation, where X]^m (k) (xL L (k')) the distorted wave in the momentum space with 

the incident (outgoing) momentum k^ (ky) and the spin projection m^. (m^^) in the asymptotic region. An important 
difference from the previous section is that the momenta k and k' are now the integration variables, whereas in PWIA 
they are fixed as k = k^ and k' = k/. To clarify this difference we use the notations 

q- = k' - k, (4.2) 

As in the previous section, we again adopt the optimal factorization approximation and set q* ~ Qc, Ripi^ — 
-^spin = 1 tt^'^s identify the [q*,n*,p*] frame with the [qc,nc,Pc] one. Then the NN ^-matrix in the optimal 
frame is written as 

q{k',k) = ^ tl,,4k',k)aoa^akb^ = J„(k',k) tl^b'{>^',>^)'TOa''Tkb' 
= JriO^' [A'cToQakO + -B'cron-fTfcn' + C"(cron* CTfco +Crooa'fc„*) 

+ E'aoq'CTkq' + F' aop' cTkp'] , (4.4) 

where a* ,b* — 0,q* ,n* ,p* and the momenta k and k' are determined by k and k' through the coordinate transfor- 
mation between the optimal frame and the NN cm. frame. 

B. DWIA t-matrix in coordinate representation 

Since the distorted waves Xkm usually calculated in the coordinate space, we now move to the coordinate 
representation. 

If t^{k' , k) were a function of only q* , i.e. f''(k', k) — t'^{q*), its coordinate representation would be local and 
would be written as 



K(k/,k.)]„^ = (xL^l. (r)*„lEv^-(r-rfc)l*noX^i.^(r)) (4.5) 

^ ' fc=i 



with 



and 



Vfe(r-r,) = I ^i^(q*)e'<i*-(-r.) (4.6) 



How ever , this is not the case in general. Neither the amplitudes t^.^. (k',k) nor the directions a*,b* — n*,p* in 



Eq. (4.4) are not determined only by q* 

The amplitudes t'^,{K',K) depend on both qc = — k = q* and Qc = /«' + k as well as the incident energy. 
Taking the same approximation as Love and Franey p5| , we suppress the Qc dependence and regard the amplitudes 
as only functions of q*, namely i™^/(q*)- We also replace the Moller factor J^(k',k) by Jrj = Jr]{kf,ki) and thus it 
becomes independent of k' and k. This approximation works for small angle scattering at high incident energy. 

As for the spin parts, the relations B' ~ F' and C" = were assumed in Ref. |l^. Then they became only q*- 
dependent, because (TQn*<^kn-' + '^ap'-'^kp* — (""o x q*) • (crj. x q*). However, the approximation becomes poorer as the 



incident energy increases |5l|,|5^. Therefore we rewrite and approximate i/:(k', k) of Eq. (^) as 
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i^(k',k)«t_^(q*) 

= [A'{q*)aooCTko + E'{q*)((To • q*)(crfc • q*) + F'{q*)((To x q*)(crfe x q*) 

+ {B'{q*) - F' {q*))aon(Jkn + C"(q*)(aonCrfco + (Toofffen)] (4-8) 



where n* in the last two terms is replaced with the averaged normal vector n of Eq. (2.5), which is independent of k 

and k'. 

We now obtain the local interaction V/c(r — r^) in the coordinate space from Eq. ( |4.6D . Consequently the matrix 
elements of the DWIA i-matrix can be written as 

[r°^(k^^k,)]„^ „^ =(*„|5 .J*„J (4.9) 

with 

5„,^,„,^ ^^(xL;L/r)|V..(r-rfe)|x^i^^(r)). (4.10) 

k 



C. Cross Sections and Spin Observables 



Inserting the DWIA i-matrix of Eq. (|4.9| ) into Eqs. (2^) and ( 2.12 ), we obtain the DWIA formulas for the cross 
sections and the sp in observables. In this procedure we need to prepare the response functions for the operator 



m^^m,^ of Eq. (^) 



2Ja + 1 



-Tr' 



n no ^ f 



(4.11) 



Explicitly writing the operation Tr in Eqs. (2.9) and (2.12), we express the cross section and the spin observables 
by Rs as 



K 



IPy = Kim [[Rs 
I Ay = Kim ([Rs 

ID,, 
ID 



-;+^ 



J++; — h 



yy 



ID, 



ID, 



+ [Rs]+^..— 

KRe{[Rs]++,__ + 
/^Re([i?5]+_.++- 
KRe{[Rs]++,__-[Rs]+^,^+ 
KRe([Rs]++,_+-[Rs]^_.^+_ 

[Rs]+-.+- - 



[Rs 



-+;-+ 



[Rs 



K (, , 



[Rs 



-+;-+ 



(4.12) 



where the suffices + and — means = ^ and — ^ respectively. 



V. ANGULAR MOMENTUM REPRESENTATION 



Let us now move on to the angular momentum representation in which calculation is actually carried out. We 
adopt the [x, y, z] frame and use the spherical tensor form cr|.^ for the spin operators and ri^ for the normal vector n, 
as follows: 
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^kO = Ifc '^k ±1 = T ^^^^^ , (^kO = ^kz 



n±i = T 



fix ± iny 



= , no = = 0. 



\/2 \/2' 

According to Eq. (fl.^), the interaction Vfe(r — rfc) can be decomposed as 

Vfe(r - rfc) = V,^(r - r^) + Vf (r - r,.) + Vf (r - r,.) 

+V,^-^(r - r,.) + V'^^ (r - rfc) + ^ (r - rfc). 

Each term of the r.h.s. is given in the angular momentum representation as 
A-term 



Vfc^(r-rfc) = J, 



^A'(q)aooafcoe''i('^-"-'=) 



{2n) 



J.-E / n{qrk)^{q)M{qr)q^dq[,%{v,)y.<jl]^^ [i'l^(f) x a 



; t 



Ey,^(rfc,r)[i'r,(ffc) X al]^ [^'Yi{v) x 



Im 



£^-term 



Vfc^(r-rfc) = J, 



^i?'(q)(<To-q)((Tfc.q)e''i('--'^) 



(27r) 



— 2 /■°° 

■/jf-X^X!"'^'"""' / 3M'rk)E'(q)jv{qr)q^dq 



M 



x[i'rz(ffc)xai]^^ [i'rp(f)xai 



JM II' 



1 J 
M ' 



F-term 



Vfc^(r - rfc) = J, j -0^F'{q){rT, x q)(^fc x q)e"i-(-'-.) 



'^'':;;:EE(^"' ~ / iiiVk)F'{q)jr{qr)q^dq 



M 



J 

A/ 



(B - F)-term 



Vf-^(r-rfc) = J„ 



1 f d'ci 



(27r) 



(B'(g)-^^'(g))ao„afc.„. 



jiq-(r-rfc) 



^ 2 



J^-EE E E^^'" 1^1 JM)(/m Ifi'l J' M')n^hl, 

^ ;m JA/ J'A/' fj.fj.' 

POO 

X / jK9^fc)(S'(g)-F'(g))jK90'z'd(7 



E E Y.yuM'rM'M[^yiiik) X m-) X -^]m^ 



x[i'll(ffc)xai],V [i'l1(f)xai];;;, 



1 J' 



JM J'M' I 
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Ci-term 

Vf^(r - r,) = ^C'(q)ao„afcoe''i ('--'=) 



/m JM 

oo 

X / ii{qrk)C'{q)ji{qr)q^dq 
x[i'y,(ffe)xa^]l'[i'rKf)xai],", 



EE^WM(^fe,^)[i'>l(ffc) X a°]l^ [i'll(f) X al\[^, (5.8) 



Co-term 



Vf^(r - rfc) = J, 1 A^C'(g)aooa,„e"i ('-"-'') 



= •^»?-EE'^(J^^-™)^'™ M -m\JM) 



TT 

/m JJ\/ 



X / ii{qrk)C'{q)ji{qr)q^dq 



= EE^™./M(^''''')[i'^'(^-^) X ^1]m [i'^'(r) X ^o"]l , (5.9) 

Im JM 



where ji{x) is the spherical Bessel function of the order I, {lmsfi\JM) is the Clebsch-Gordan coefficient and aji = 
(J 1 0|^ 0). The spherical tensor product is expressed as [i'Y; x cr"]^^ = X]m^i(^'^*Ml^-^^)i'^m(^)o'^- 

For an optical potential with a spin-orbit force, the distorted wave Xu-ii (i") is expressed in the angular momentum 
representation as |Q 

= "T" E E (^^^i^^sMimj.)yCmS^i) 

x(Z,m^sm;jj,m,jy,,,„,(f)z''e''^'.Mg(fc,,r)|m;^), (5.10) 



where \ms) is the intrinsic spin state with the z-projection rris and ai is the Coulomb phase shift. The radial part of 

4/ 



the distorted wave u[^\k,r) has the asymptotic behavior 



u'if{k,r) ~ e**'^ sin(fcr- 7/cln2fcr- ^-^ + cri + Sij), (5.11) 
■' 2 



where rjc is th e So mmer feld p arameter and dij is the nucl ear phase shift. 
Using Eqs. (^.3D and (5.10) together with Eqs. (^.4D - (5.9), we can write as 



where 



(xK,^.(r)|V,(r-r,)|x^ijr)) ^ ^ 5(7;/,";* (k,, k,; r,) [i^yK?.) x <jIYJ , (5.12) 

IsJM 



^(JO)JM (k/, k,; rk) ^ r drVj {rk,r)f^j^'^jj^; (k/, k,; r) 
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/•OO 

Im ^ 
nOQ 

/•OO 

^ "'0 



with 



/247r 



xu|+^(fc/,r)w^(fc^,r)(2^, + l)y^(2ji + l)(2;/ + l)(Z,0Z0|Z/0) 
s/ J/ 

X Si ji ^ (?i0sjmsJjimsJ(//r7i/S/TOs^|j7ms, + M) 

Z s J 



We can now write Sm^^m,. of Eq. ( 4.10| ) as 



where 

A 



PiuUMir) = E ^^^T^if^ l^'Y.ih) X <5 f ^ r,, ) , (5.18) 

A:— 1 Vfc'^l / 

which is the radial part of the transition density 

k=l \k' = l / 

= {lmstx\JM)p(,,)jM[r) [I'YUv)]* . (5.19) 

IsJM 

VI. CALCULATION OF RESPONSE FUNCTIONS 



(5.13) 



(5.14) 



/(T;)SM"(k.,k/;r) (5.15) 



x(j-m,.JM|j/m,, +M)(-)^^^^y'|^^-^h^ (5.16) 



EnOO 
/ '5(;.)jj\/'(k»,k/;r)pJ;^)^j,,^(r)r dr, (5.17) 



In Sec. IV, we presented a method of calculating the cross sections and the spin observables by using DWIA. It 
was shown that what we need are the response functions Rs{lu) of Eq. ( 4.11| ). Here we describe a way of calculating 
them. 
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A. Polarization propagator 



Following Refs. |Q and [g^, we introduce the polarization propagators for the spin-dependent transition density 
in the angular momentum representation as 



'^{i's')j'M';(is)jM{'r' , r; co) 



1 



2Ja 



n no 



\p(i's')j'M' {r')\ *„> I Pli,) jM (r) I ) 



IT] 



(*"olP(,,),/MWI*")(*n|P(('.')^'A/'(^')l*no> 



2Ja 



-LU - {£n - £no) + 

P{i's')j'M'{r') 
1 



LU - {Ha ~ £no) + "7 

l*no) 



p\ls)JM^''') 



+Piis)jMir) _^ _ _ ^^^^ ^ .^P^i's'U'M'ir') 



(6.1) 



By using the Wigner-Eckart theorem, we can prove that ^{iis')J'M';{Is)jm £^re diagonal with respect to J and M and 
independent of M. The corresponding response functions are defined as 



Rj(l's')ils){r',r;uj) = ^ ^'(*„olP(/'.')./M('^')l*n>(*„|p{,,)jM WI*"o>^ - (^n - ^"o)) 

n no 

^ -■^^^^J{i's'){is){r ,r;uj). (6.2) 
Using Eqs. ( p7| ) and ( ^l8| ), we obtain the response function Rs{oj) of Eq. ( [4.11| ) as 

poc poo I I 



^Rj{i's')iis) {r', r; w)5(™^j7/' (k/, k^; r). 



(6.3) 



Inserting these Rs into Eq. ( 4.12| ), we can calculate the double differential cross section and the spin obsevables. 

In the following we consider only the case where the target nucleus is a doubly closed shell. Therefore we can set 
Ja = and omit the summation over rig, and denote £„(, and simply by £o and , respectively. 



B. Single particle model 

A simple approximation to evaluate the polarization propagators is the single-particle model. In this model Ha is 
replaced by the single-particle-model Hamiltonian 

H^'^=T.(-^ + U')+^Eo (6.4) 

where to — [(^ — l)/^]TOAr, Uk is the mean field for the /c-th nucleon, and Ai?o consists of the masses of the constituent 
nucleons and the energy shift to give the correct total Hartree-Fock energy of the target ground state. Since there 
appears only the energy difference between the target ground state and the final state of the residual nucleus for the 
response functions, we discard the term A_Bo from now on. 

The mean field is determined as the Hartree-Fock field for the target nucleus and the single particle energy 
is given by the equation 

f-J^+f/^a =ea0a. (6.5) 

We write Cq, of the unoccupied state by Cp and that of the occupied state by e^. Note that is measured from 
one-nucleon separation threshold of the target. 
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1. Polarization propagator 



The polarization propagator is given by 



n 



ph 



(^o|p(i'.')JM('^0bfe)(pfe|4,)jM(Ol^o) 



+ 



14.) JM (^) \P^) WP{1's')jm {r) I $0) 



-uj - (Cp - eh) + IT] 



(6.6) 



where \ph) denotes a one-particle-one-hole state and $0 is the ground state of the target in the single particle model. 
This is called the uncorrelated polarization propagator. 

The infinite sum over p can be handled by the single-particle Green's function as 



ph 



B[vs')j{h,p)— — —B^i^-^j[h,p) 



+ I3{1'S')J[P, — — h) 



(6.7) 



where p and h denote the sets of the single-particle quantum numbers p — (IpSpjp) and h = {nhlhShjh), respectively, 
and 



ja 



S(z«)j(a,6) = V2ja + lV2jb + l<( lb Sb jb } {la\\i'Yi\\h){sa\\a'\\sb). 

I s J 



(6.8) 



The radial part of the single-particle Green's function is given by 

2m 



gp{r',r;E) 



Wifp,hp) 



Jp{r<,E)hp{ry,E) 



(6.9) 



where W is the Wronskian and ?'<(>) denotes the smaller (larger) one of r an d r'. The radial wave functions Uh, fp 
and hp are the bound state, the regular and the irregular solution of Eq. ( |6.5|) , respectively. 



2. Spreading widths 

To take account of the spreading width of the particle states, we adopt a complex potential for the single-particle 
potential U as in Refs. 



U{r) = V{r)+iW{r), 



(6.10) 



while we use the real potential for the hole states. For such a choice, the orthogonality between the hole and the 
particle wave functions is destroyed. We therefore utilize the orthogonality condition prescription |l^ , details of 
which are explained in Ref. |Q . 

To take account of the spreading width jh of the hole states, we replace the hole energy eh by the complex energy 



eh = e/i + 1- 



.7h 



(6.11) 



3. Effective mass 



In principle the single-particle potential U is non-local, and thus the Schrodinger equation for the single-particle 
states is 
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- 7r-V^<?!>(r) + / dr'[/(r,r')(/)(r') = e0(r). 
2m J 



(6.12) 



We deal with this non-locahty by an effective mass approximation . We first introduce the Wigner transforma- 

tion of the potential U{r,r') as 



[/(r,q)= / dsC/(r+|,r-|)e''^•^ 



(6.13) 



We assume that U is spherically symmetric with respect to r and q, respectively, and write it as U{r, q^). We expand 
it with respect to around the square of the local momentum K(r) 

n^{r) ^2me-U{r,K^{r)), (6.14) 

and we keep only the terms up to the first order. Introducing the radial-dependent effective mass m* (r) as 



!^^ = /(.) Jl + 2m^%^ 
m \ OK'' 



(6.15) 



we obtain the equation 



2m* (r) 



where the local potential Ui^{r) is 



ULir) 



1 



8 V ■m.*{r) 



(r) = e(j>{r), 



Ui^{r) = —{U{r,K\r))~e) + e. 



(6.16) 



(6.17) 



Since we do not know about U{r,r'), we determine Ui,{r) and /(r) phcnomcnologically. Note that here we only 
consider the fc-mass but not the E'-mass. 



C. Ring approximation 

Next let us take account of the nuclear correlation by means of RPA. To treat the exchange terms of RPA we 
further utilize the ring approximation, which replaces the exchange effects by a contact interaction. The polarization 
propagator then satisfies the RPA equation p3| , p6| , p7| 

+ ^ '"i^'"i''2f^''2n^°Lo(hsi)(''''''i;'^)^-f('i^i)('2S2)(''i'^2)nj(,3^3)(;^)(r2,r;w) (6.18) 



hsil2S2 

where VFj(;jsj)(;2S2)(^i: ''2) is the radial part of the effective interaction as 

V{rur2)=J2 E [i'^>^;.(fi)xaJ^]i;W^j(i,.,)(i,.,)(ri,r2)[i'=rz,(f2)xa2^^]i,. (6.19) 

JAI I1S1I2S2 



D. Isospin and A degrees of freedom 



From Sec. IV up to here we have suppressed the isospin and the A degrees of freedom, but of course we have to 
include them in the a ctual calculation. Here we summarize the final formulas including these degrees of freedom. As 
was discussed in Sec. HID, the relevant quantities carry the superscript a {— A^, A, A*) and the isospin subscript n. 



The extension is straightforward a nd ea sily understood. 
The transition density operator ( 5.18| ) is to be 
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Pk(Is)JM 



k=l 



rrk 



(6.20) 



and thus the response function (3.2) becomes 



n no 

X,5 {UJ - {£n - £no)) ■ 



(6.21) 



Further details are given in Ref. 



The NN t-matrix was generahzed as Eq. ( 3.24 ) . Correspondingly, the source function ( 5.12 ) is generalized as 



JV(+) 



J t 



(6.22) 



a IsJM 



with N,N' — p or n, which specifies the reaction type. Here the generalized forms of Eqs. ( 5.13|) and ( 5.14 ) are 



with 



'5.(i.)SA/(^ k/,7Vk,;rfe)-^ / r'drV:.''{rk,r)f... ' ' (TV'k/, TVk,; r) 



rns ^ ras , , 

/(z.)SM'(^k^^iVk,;r) 

= E / dl^rX™^i;.(k,,r)«^| [z'y,(f) X or\{, |mUx™[1..^(k.,r), 

where V^.^ are the generahzed forms of Vf, Vi^j^, Vfii,,Vfii, and Vi%^j,^.j,. 
Finally, we obtain 

L J m= , TUs ■ ;m' m' 

Sf s^: Si 

^(iV|r„7|A^')(iV'|ro^J^>EEEE/ '"'^'\ ^''^r' 

8m' m' * a am^ m« 

S^v^jM {N%,N\.,- r')<,(,,,)(,,) (r', r; ^)5^(,,)S^' (iV'k^, TVk,; r). 

in place of Eq. (|6.3|) . 



(6.23) 



(6.24) 



(6.25) 



VII. NUMERICAL CALCULATIONS 



We apply the DWIA formalism presented in the previous sections to the intermediate energy (p, n) reactions around 
the quasi-elastic peaks. We calculate the polarized cross sections IDq and I Dp of the reactions observed at LAMPF 
§ and RCNP @ as summarized in Table |. 



A. Choice of parameters 

First we present the parameters that we selected for our numerical calculations. 
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1, Optical potentials 



For the optical potentials of the incident proton and the outgoing neutron, we took the global optical potential given 
by Cooper et al. It is designed for the proton but we used it even for the neutron by excluding the Coulomb 

part. To study the optical potential dependence, we also used the proton optical potential given by Jones et al. 
p3[ and the neutron global optical potential given by Shen et al. . For the potentials which have the form of the 
Dirac phenomenology, we rewrote them in the Schrodinger equivalent form and solved the non-relativistic Schrodinger 
equation to obtain the distorted waves. We used the relativistic kinematics and the reduced energy prescription. 



2. Effective mass, single-particle potentials and spreading width. 



To obtain the effective mass, we have to determine /(r) and UL{r), as was discussed in Sec. VI B 3. Suggested by 
Ref. Esl, we assume f{r) has the form 



fir) 



•m*{r) 



1 



1 + exp I 



with R = rpAi, since we consider only the k-mass. Note that m*{r — 0) 
to reproduce the energy spectra of ID^. 
For Ul we used the form 



(7.1) 



(1 — h)m. The parameter b will be adjusted 



UL{r) 



-(K) + iWo) 
2 Vis 



l + exp(r-ii) 
exp(r^) 



mt a 



(l + exp(l^)) 



1 • s + -t^c. 



Oul ; 



(7.2) 



where Vcoui is the Coulomb potential of a unifor mly charged sphere with the radius Rc = TcA's . The shape parameters 
are fixed at = = 1.27 fm and a = 0.67 fm |47]| , and the nucleon spin orbit potential depth Vis is chosen to be 6.5 
MeV for ^'^C and 10.0 MeV for ""^Ca, respectively The real potential depth Vb is determined in such a way as to 
give the observed separation energy of the outermost occupied state of the target nucleus. 
Using the phenomenological energy-dependent relation for the spreading width Bst] 



7(e) = 21.5 



(e-ep) 



(e-eF)2 + 182 



1102 



(e-eF)2 + 1102 



(7.3) 



in units of MeV with the Fermi energy ep, we set the imaginary potential parameter Wq for the particle to be 
Wq = \l{^p)i ^iid the hole spreading width 7/1 to be 7/j = 7(eft,). 

For A, we set mA = 1232 MeV, Vb = 30 MeV and Wo = Vis= 0.0 MeV [p7|. 



3. Effective interaction for ring approximation 

We take RPA correlation into account only in the isovector spin- vector channel (T = 1, 5 = 1 in Sec. |IIID| ). 
For other channels, we simply use the uncorrelated response functions. For the effective interaction in the isovector 
spin-vector channel, we employ the (vr + p + g') model, in which it is written as 



l/"^(ri - ra;^) = (r? • r^) \v^^{vi - v^-.io) + Vf {v^ - v^^oj) 



with 



C^(ri-r2;t^) 



^^'^(ri -V2\uj) 



(27r) 
(27r)3 



(7.4) 

(7.5) 
(7.6) 



Here W^'^ and W^'^ are given by 



Q/3 
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Wl^qM = ^ {g'^, + Tl[q^u)Tl{q,u.)-^-^—-^ ) , (7.7) 



^ - 9 - ™^ 

2 



w^'{<i^-) - ^ (.9l, + c^,rs(9,c.)r^(g,c.)^^^-^-^) , (7.8) 

where m^^ {— 139 MeV) and rup (= 770 MeV) are the pion and the p-meson masses, respectively, and the coefficient 
C^p is the ratio of the p-meson coupling to the pion one. The g-dependence of the Landau-Migdal parameters g'^^ is 
neglected and the vertex form factors are chosen to be 

_ A2 7T)2 _ A2 

Vl{q,L0)= , % A2 ^ r^(g,c.)^ ^ ^ , " (7.9) 

where A^r and A^ are the cutoff parameters. We note that we can identify A* with A in 5^^, /„, C^^ and T^''\ 

We set the parameters as = 1.0, /a/Zw = 2.0, C^^ = 2.18, A^ = 1300 MeV, and Ap = 2000 MeV §. As to the 
Landau-Migdal parameters, we fix g'j^^ — 0.6, = 0.5 and we adjust to reproduce the observed IDq as well 
as possible. 



^. Energy Shift 

The nuclear structure model used here is a simple one which is based on the mean field approximation with the RPA 
correlation, because we are interested in only gross structure of highly-excited states in the continuum. Apparently 
the model does not well represent the structure of the low lying states, especially the ground states of the deformed 
nucleus ^^C. Therefore the excitation energy obtained by the present model should be smaller than the observed one, 
because the real ground state energy should be lower than that given by the model. 

To remedy this shortcoming, we made the following modification. The observed energy spectra of ^^C show 
the eminent peak of the unresolved 4~ and 2~ states of ^^N. We artificially add 5 MeV to the energy transfer lo to 
coincide the calculated peak to the observed one for the uncorrelated case (to* = to and no RPA correlation). For 
the correlated case, this energy shift becomes 3.5 MeV, but we kept the 5 MeV shift for all different parameter sets. 
This small difference of energy shift does not affect our conclusion. Any shift was not made for ''"Ca. 



5. Convergence 

In the actual calculations, we have to limit the infinite range of the summations and integrations by judging from 
their convergence. The maximum angular momenta of the distorted waves were fixed at 40, and the maximum 
transferred angular momenta were set at 9 for ^^C and 12 for ^"Ca for both the incident energies, because the 
transferred momenta are close in both cases. 

The maximum momentum for the Fourier transformation from the momentum space to the coordinate space was 
set to be 4 fm~^, though the main contribution comes from q near the observed transferred momenta shown in Table |. 
The radial integration was limited up to 10 fm, since the response functions are well damped at this radius. 



B. Results of DWIA calculations 



First we try to reproduce the spin-longitudinal cross section IDq of the ^^C(p, n) reactions at 494 MeV by adjusting 
the Landau-Migdal parameter and the effective mass at the center m*(r = 0), namely parameter 6, within a 
reasonable range. 

The results are shown in the IDq part of Fig. [T]. The dashed line denotes the result with to*(0) = to (i.e. 
h — Q) and without the RPA correlation. The dotted line denotes the RPA results with the universality ansatz 
9'nn ~ 9na ~ 5aa ~ 0.6 and again with m* (0) — to. We see that the experimental data are much larger than these 
two results, and therefore we reduced g^^, guided by the fact that the response functions increase as .g^^ decreases 
pTf . The dot-dashed line shows the RPA results with (s^vTV' 5jvA' 5aa) ~ (0.6,0.3,0.5) and m*(0) — to. Now IDq 
drastically increases but the peak position is much lower than the observed one, owing to the softening. We then 
reduced the effective mass to*, guided by the prediction of the Fermi gas model (wpcak = q'^/2m*) and the numerical 
calculation of Itabashi The full line represents the RPA result with (s^vTV' SwA' Saa) ~ (0.6,0.3,0.5) and with 
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m*(0) ^ 0.7m(6 = 0.3). By this change of m*, the peak position moves up very close to the observed one, though the 
peak value is reduced. This is consistent with the Fermi gas prediction (3m*/4qpF, Pf being the Fermi momentum). 
Now the fit is very much improved by the choice of these parameters for IDq. 

The results of IDp are shown in the right panel of Fig. |^. We see that the RPA correlation with 5^jv = 9'na ~ 
g'^^ = 0.6 markedly quenches IDp, as was predicted, whereas the observed data are very much enhanced. When we 
reduce g'j^^ the results increase considerably, but are still quenched at low lo. When we further reduce the effective 
mass, the peak shifts upwards. We found that all of the calculated results are much smaller than the observed IDp. 
In the end, we could not reproduce the observed IDp by changing these parameters within reasonable ranges. 

With the same sets of parameters we calculated the other reactions given in Table |. Fig. |^ shows the results for 
'^°Ca{p, n) at 494 MeV and Fig. | shows the results for ^^C, ^oca (p, n) at 346 MeV. The features of the results are 
common for all cases, though the fit of IDq is somewhat better for the case of 494 MeV than for that of 346 MeV. 



VIII. TEST OF EFFECTIVE NUCLEON NUMBER APPROXIMATION 



In the previous experimental papers 1|M| the response functions were extracted by means of P WIA with the effective 



nucleon number approximation described in Sec. III. We test this prescription by comparing the results with those 
of DWIA. 

For detailed analysis, we also introduce mode-dependent effective nucleon numbers defined by [|j 

A^c'-ffH-A^ ^PWj^PW(^) , ^cy-) = ^7fnv^|w|^- (8-1) 

They depend on the spin-longitudinal (L) and spin-tra nsver se (T) modes as well as the excitation energy w, though 
the effective nucleon number A^ff determined by Eq. ( 3.1^ ) (effective neutron number A'eff in the present case) is 



independent of them. 

In the present analysis, we try to get N^^ [uj) without the ijj-dependence by estimating them at a certain to near 

the peak. The obtained n'^^^^ with and without the RPA correlation are summarized in Table ^ 

For the uncorrelated cases, we found that iV^ « A^^g and thus we set them equal and denote both of them by 

n'^ in common. We note that Wakasa et al. |^ reported that A^ cff es timated by Eq. ( [3.1^ ) and n'^ obtained above 
are very close, and also note that the estimation of N^s by Eq. ( p.lg ) is somewhat ambiguous because the equation 



includes the uncertain total A^A^ cross section in the nuclear medium. Therefore here we identify A^eff with N^^ . 

The PWIA cross sections multiplied by N^ajN — Ncff/N (dashed lines) and the DWIA cross sections (solid lines) 
without the correlation are compared in Fig. ^(a) and (c) and Fig. ^(a) and (c) for cases of ^^C and ""^Ca with the 
incident energy 494 MeV, and ^^C and ''"Ca with 346 MeV, respectively. The ^-independent approximation holds 
well only for the IDq of ^^C at 346MeV, but not so well for other cases. This is an indication of poorness of the 
effective nucleon number approximation. 

The results with the RPA correlation are compared in Fig. |J(b) and (d) and Fig. ||(b) and (d) for cases of ^^C and 
^°Ca with 494 MeV, and ^^C and '^°Ca with 346 MeV, respectively. The DWIA results, the PWIA results multiplied 
by N^g ' /N and the PWIA results multiplied by Nes/N are denoted by the soHd, the dashed and the dotted hues, 
respectively. The RPA calculations are carried out with (ff^vAr; 5jva> 3aa) ~ (0.6,0.3,0.5) and to*(0) — 0.7m. Once 
the RPA correlation is included, Njf^ is increased but A^ is reduced and thus they differ very much to each other. 
Especially the difference is very large for ''"Ca. This is another strong indication of poorness of the Acff method. The 
above changes of N^g and A^ due to the RPA correlation may be explained in the following way. The enhancement 
of i?L and the quenching of Rt are stronger in the higher-density region (the inner part) but they are masked by the 
stronger absorption in this region. Consequently, we should decrease A^ to reduce the enhancement seen in PWIA 
but increase Nj^ to reduce the quenching. 

The large difference, especially for ^''Ca, between the results of the effective nucleon number approximation (the 
dotted lines) and the DWIA results (the solid lines) clearly shows that the conventional way of obtaining the response 
functions does not work quantitatively. Therefore the DWIA calculation is definitely necessary in the quantitative 
analysis. 



IX. DISCUSSION 
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1. Effect of the spin- orbit force 



One of our main concerns in the DWIA calculation is the effects of the spin orbit force of the optical potentials. 
We compared the numerical results with and without the spin orbit force in Fig. ^ for ^^C(p, n) at 346 MeV. The 
RPA correlation was not included in this analysis. We found that effects are larger for D^"^ than for D^"^ . 
Fortunately, however, they are so small that the spin-orbit force does not greatly disturb the separation between the 
spin-longitudinal and spin-transverse responses. We also reported similar results for ^^C(p, n) at 494 MeV p9| . 

2. Ambiguity of the optical potential 

To investigate the effects of the ambig uity of the optical potentials, we compared the results in terms of the global 
optical potential given by Cooper et al. ||42f| with those in terms of the proton optical potential given by Jones et al. 
p3[ and the neutron global optical potential given by Shen et al. Q]. Fig. |^ shows the comparison of Ii^"^ and 
/°^Dp ^ for ^^C(p, n) at 346 MeV. Here the RPA correlation was not included. We found that the ambiguity hardly 
affects the results. 

3. Ambiguity of the NN amplitude 

It has been pointed out Q that the theoretical results are affected to some extent by the ambiguity of the NN 
scattering amplitudes obtained by the phase shift analyses. To determine the influence of this ambiguity we compared 
the results with the NN amplitude obtained by Bugg and Wilkin Bl| and with that obtained by Arndt and Roper 
||. Fig. I shows the results of I^^^ D^^ and I^^^ D^^ for ^^C{p,n)at 346 MeV without the RPA correlation. The 
effects are less than 10%, and they are different for /^^IJ^^ and I^^^ D^"^ . 

4- Two-step processes 

De Pace calculated the one- and the two-step contributions in terms of the Glauber approximation. He found 
that the two-step process is more effective for I Dp than for IDq. He concluded that IDq is reasonably explained 
but that there is a large discrepancy between the observed and the one-step results for I Dp, and that the two-step 
contribution is sizable but not sufficient to provide an explanation of the large discrepancy. The A contribution was 
not included in his analysis. He made the simple assumption for the two-step process that the first step is caused by 
the spin-scalar isospin-scalar and the second step by the pure E (F) amplitude and vice versa. 

Nakaoka and Ichimura ||3l|| evaluated the two-step contribution through its ratio to the one-step one obtained by 
PWIA with an on-energy-shell approximation. They took into account all five terms of the NN amplitudes in both 
the first and second steps. They also found that the two-step process is more effective for I Dp than for IDq. For 
IDq it accounts in large part for the discrepancy between the DWIA and the experimental results in high energy 
transfer regions. They also found that the two-step contribution does not help to reduce the large discrepancy seen 
in IDp. Recently, Nakaoka reported that the two-step effect becomes twice as much as the previous result, when 
the off-energy-shell effect is included. 

5. Remaining problems 

There still remain various other points to be considered, especially keeping in mind the large discrepancy in IDp. 
On the structure side, they are (1) nuclear correlations beyond RPA which has benn evaluated by various authors 
[ p^ , ^JT2| , (2) use of genuine RPA Q instead of the ring approximation, and (3) self-consistency between the mean 
field and RPA | ]57| , |58[ , etc. On the reaction side, they include (1) removal of the optimal factorization approximation, 
(2) removal of the averaged fi approximation, (3) off-shell effects of the NN ^-matrix, and (4) more realistic 
NN — NA transition t-matrix, etc. 
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X. SUMMARY AND CONCLUSION 



We first summarized the general formulas for the polarized cross sections, IDi, of the {N,N') scattering to the 
continuum, which are closely related to the nuclear spin response functions. We then reviewed the PWIA formalism for 
the intermediate energy (iV, N') reactions, treating the Fermi motion by the optimal factorization and the absorption 
by the effective nucleon number approximation. The A isobar degrees of freedom are also included. This has been 
widely used to extract the nuclear response functions. 

Since the reliability of this method has been questioned, we here developed the DWIA formalism for the intermediate 
energy (iV, N') reactions to the continuum incorporated with the continuum RPA method. This formalism is described 
in the angular momentum representation in the coordinate space with the response functions, Rj(v s'){is){'i'' 
which are non-diagonal in the coordinate space as well as in the spin space. The formalism still involves the optimal 
factorization and the averaged reaction normal approximation. 

The response function calculation includes the radial-dependent effective mass and the spreading widths of the 
particle and the hole. The RPA correlation are calculated by use of the interaction of the {n + p + g') model. 

The presented formalism is applied to ^^C, "^^Ca {p,n) at 346 and 494 MeV in the quasi-elastic region. In this 
analysis, the Landau-Migdal parameters and the effective mass are treated as adjustable parameters. 

From this analysis, we draw the following conclusions: 

(1) The observed spin-longitudinal cross sections IDg are reasonably well reproduced by adjusting g'j^^ and m* and 
by adding the 2-step contribution. The analysis indicates that the smaller g'j^/^ (~ 0.3) and the smaller effective mass 
at the center (m* « 0.7) are preferable. This claim of the smaller is consistent with the conclusion obtained from 
the analysis of the GT resonances by Suzuki and Sakai . 

(2) The observed spin-transverse cross sections I Dp are much larger than the calculated ones and cannot be reproduced 
within the present theoretical framework. However, one must note that Rt corresponding to the observed IDp is also 
much larger than that obtained by (e,e'). This confirms the previously reported conclusion [^J^. 

(3) The iVcff method, which has been conventionally used to obtain the response functions, is found to be a quanti- 
tatively poor approximation, especially for heavier nuclei such as Ca. The DWIA analysis is definitely needed in the 
quantitative analysis. 

(4) From these findings, we stress that the observation oi Rl/Rt < 1 does not necessarily constitute evidence against 
the enhancement of Rl and of pions in nuclei. We rather conclude that the present experimental data support the 
enhancement of Rl- The contradiction regarding the ratios comes from the extraordinarily enhanced IDp. Before 
drawing any definite conclusion about the response functions, we must disentangle the anomalously large IDp puzzle. 

(5) We also investigated the effects of the spin-orbit distortion and the ambiguity of the optical potential and the NN 
i- matrix, and found that they are not especially significant. 
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TABLE I. Reactions to be analyzed. The values of gem are those around the peak. 



24 



Target g..iab (MeV) ^lab gem (fm"^) Refk 





494 
346 


18° 
22° 


1.70 1 
1.67 








494 
346 


18° 
22° 


1.71 1 
1.68 1 







TABLE II. Effective neutron number. "No" means the no-correlation (uncorrelated) cases. 



E, (MeV) 


Target 


Correlation 






494 




No 


2.17 


2.17 






RPA 


2.12 


2.39 




«Ca 


No 


4.00 


4.00 






RPA 


3.60 


5.78 


346 


V2^ 


No 


2.58 


2.58 






RPA 


2.53 


2.82 




*OCa 


No 


4.62 


4.62 






RPA 


4.25 


7.00 
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FIG. 1. DWIA results of IDq and IDp for ^^C(p,n) at 494 MeV. The dashed hne denotes the result with m* = m and 
without the RPA correlation. The dotted and the dot-dashed lines represent the RPA results of (<?jvjvi 5jvai 5aa) = (0-6, 0.6, 0.6) 
and (0.6,0.3,0.5) with m* — m, respectively. The solid line shows the RPA results of (0.6,0.3,0.5) with m*(r = 0) = 0.7m. 
Experimental data are taken from Taddeucci et al.[8] 
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FIG. 2. The same as Fig.0,but for *"Ca(p,7 
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FIG. 3. The same as Fig.§, but for (a) "C(p,n) and (b) '^°Ca{p, n) at 346 MeV. The experimental data are taken from 
Wakasa et a/. [91 
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FIG. 4. Comparison of DWIA and PWIA results at 494 MeV: (a) without and (b) with the RPA correlation for ^^C(p,n) 
and (c) without and (d) with the RPA correlation for *''Ca(p, n). The solid and the dashed lines denote the DWIA and the 
PWIA results multiplied by Nj,^^ /N{^ Ncs/N) in (a) and (c) and by N^'^'^^/N in (b) and (d) , respectively. The dotted lines 
denote the PWIA results multiplied by Ncb/N in (b) and (d). The values of N^g and N^g are shown in the figure. 
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FIG. 5. The same as Fig.W, but at 346 MeV. 
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FIG. 6. Comparison of the DWIA results with (full line) and without (dashed line) the spin orbit force for ^'^C{p,n) at 346 
MeV. 
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FIG. 7. Comparison of the results with various optical potentials for ^^C{p,n) at 346 MeV. The full lines denote the results 
with the potential obtained by Cooper et al. [42] for both the proton and the neutron, while the dashed lines denote the results 
with the potential by Jones et a/. [43] for the proton and by Shen et a/. [44] for the neutron. 
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FIG. 8. Comparison of the results with the NN amplitudes given by Bugg and Wilkin (full line) [51] and by Arndt (dashed 
line) [52] for ^^C(p,n) at 346 MeV. 
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